Specific Heat Study of the Magnetic Superconductor HoNi 2 B 2 C 
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The complex magnetic transitions and superconductivity of H0N12B2C were studied via the de- 
pendence of the heat capacity on temperature and in-plane field angle. We provide an extended, 
comprehensive magnetic phase diagram for B \\ [100] and B | [110] based on the thermodynamic 
measurements. Three magnetic transitions and the superconducting transition were clearly ob- 
served. The 5.2 K transition (Tjv) shows a hysteresis with temperature, indicating the first order 
nature of the transition at B = T. The 6 K transition (Tm), namely the onset of the long-range 
ordering, displays a dramatic in-plane anisotropy: Tm increases with increasing magnetic field for 
B || [100] while it decreases with increasing field for B \\ [110]. The anomalous anisotropy in Tm 
indicates that the transition is related to the a-axis spiral structure. The 5.5 K transition (T*) 
shows similar behavior to the 5.2 K transition, i.e., a small in-plane anisotropy and scaling with 
Ising model. This last transition is ascribed to the change from a* dominant phase to c* dominant 
phase. 



I. INTRODUCTION 

An understanding of the interplay between mag- 
netism and superconductivity has been an area of in- 
tensive research because of their seemingly antagonistic 
tendencies J»£ The limited availability of suitable exam- 
ples and the low magnetic transition temperatures of 
those that exist have made the studies very difficult ex- 
perimentally The recently found RNI2B2C family where 
R is rare-earth element offered a new venue because the 
magnetic transitions (Tm) occur in an easily accessible 
temperature range with a variation of Tm/Tq ranging 
from 1.75 for Dy to 0.14 for Tm. Further, high-quality 
single crystals have become available. Among the mag- 
netic members, H0M2B2C is particularly interesting be- 
cause its complex magnetic phases are observed to coex- 
ist with superconductivity. Neutron scattering revealed 
three types of magnetic order A commensurate anti- 
ferromagnetic structure is formed with q = c* , i.e. (0 1) 
below 6 K in which the spins order ferromagnetically 
within the a — b plane and antiferromagnetically along 
the c— axis. Two incommensurate structures coexist over 
a finite temperature range, i.e., 5 K < T < 6 K: the 
c* structure with q — (0 0.915) where the ferromag- 
netically aligned spins are rotated by 165° and the a* 
structure with q = (0.585 0) for which the detailed 
structure is still unknown. 

H0M2B2C exhibits a near-reentrant superconductiv- 
ity, i.e., reentrant resistive behavior in a small magnetic 
field, due to competition between superconductivity and 
exchange-coupled antiferromagnetic order & In a theoret- 
ical analysis of the interplay, the onset of the c— axis in- 
commensurate state has been shown to suppress super- 
conductivity, leading to the near-reentrant behavior. 7,8 



Kreyssig et al., however, found that in Y-doped quasi- 
quartenary compound, Hoi_ x Y :c Ni2B2C, only the in- 
commensurate a— axis feature remains in the same tem- 
perature range as the near-reentrant behavior while the 
c— axis spiral structure exists over a much wider temper- 
ature range, indicating that it is the a*-structure con- 
nected with Fermi surface nesting that enhances the pair- 
breaking effect £ For a further understanding of the inter- 
play, it is necessary to study the nature of the magnetic 
transitions and to establish the B — T phase diagram of 
HoNi 2 B 2 C. 

The magnetic phase diagram of HoNi2B2C has been 
studied extensively in the context of the interplay be- 
tween magnetism and superconductivity, mostly below 
B = 1 TiS Zero-field specific hea1»i£ and other surface 
sensitive measurementsiii^ indicate that there are three 
magnetic transitions (Tjv, T* , Tm), while other specific 
heat dat a 13 i 14 show only two transitions (Tjv, Tm), where 
Tm = 5.2 K is the Neel temperature, Tm — 6 K is the 
onset of a long-range magnetic ordering, and T* = 5.5 K 
is ascribed to some change of the oscillatory magnetic 
state. It has been suggested that the low temperature 
physical properties of polycrystalline H0M2B2C depend 
on thermal treatment as well as chemical composition 
while the properties of single crystals are relatively less 
affected^£*i£ Even though majority of measurements in- 
dicate the presence of T* at 5.5 K, thermodynamic mea- 
surements on single crystals that show the magnetic tran- 
sition to be a truly bulk property are rare. 

In this paper, we report specific heat measurements 
of single crystal H0M2B2C as a function of temperature 
and magnetic field to provide an extended (up to 6 T), 
comprehensive magnetic phase diagram for B || [100] and 
B || [110]. Three magnetic transitions are clearly re- 
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FIG. 1: Top panel: the specific heat of sample A in 
zero field. Bottom panel: the specific heat discontinuity, 
AC(B,T) = C(B,T) - C{0.2T,T), in the vicinity of T c for 
5 = 0, 50, 100 mT 



solved, confirming that T* is a bulk magnetic transition. 
The 5.2 K transition (T/v) shows hysteresis with tempera- 
ture, indicating a first order antiferromagnetic transition. 
The 6 K transition (Tm), namely the onset of long-range 
order, displays a dramatic in-plane anisotropy: Tm in- 
creases with increasing magnetic field for B | [100] while 
it decreases with increasing field for B || [110]. The 
anomalous anisotropy in Tm indicates that the transi- 
tion is related to the a-axis spiral structure. The 5.5 K 
transition (T*) shows a small in-plane anisotropy and 
was ascribed to the transition from a*-dominant phase 
to c* dominant phase. The specific heat jump related to 
the superconducting transition was observed at 8 K. 



II. SUPERCONDUCTIVITY 

We have used two single crystals of H0M2B2C, labeled 
sample A and sample B for the heat capacity measure- 
ments. Both samples were from a same batch. A high 
temperature flux method was used to grow the sam- 
ple using NI2B as a solvent. The details are described 
elsewhere 4i 

The heat capacity measured by ac calorimetrjsi^ was 
converted to absolute values using literature data^ The 
upper panel of Fig. 1 shows the specific heat of H0NI2B2C 
vs. temperature at zero field. The three sharp peaks 
mark the magnetic transitions noted above. The first 



transition at 6 K is marked as Tm, the second at 5.5 K, as 
T*, and the third at, 5.2 K as T/y. There is also a specific 
heat discontinuity around 8 K related to superconducting 
transition but it is too small to be seen on the same scale. 

The bottom panel shows the specific heat discontinuity 
related to the superconducting transition in the vicinity 
of T c . The circles describe the specific heat difference, 
AC = C-C n , for B = 0T, the triangles for B = 50 mT, 
and the crosses for B = 100 mT, where C n is the specific 
heat in normal state. Data taken at 0.2 T were used as 
the normal state background between 7 K and 9 K be- 
cause the magnetic contribution is unaffected at low fields 
and the superconductivity is suppressed in that temper- 
ature range. In zero field, the transition from normal to 
superconducting state occurs at 8.04 K with a narrow 
transition width (AT/T C < 0.04). The transition tem- 
perature was defined as the midpoint of the transition, 
which is essentially equal to that found from entropy- 
conserving method. As the magnetic field increases, the 
T c decreases at the rate -8.3 K/T and the transition width 
becomes broadened. 

The specific heat jump in zero field is about 110 
mJ/mol-K. If we use the BCS relation AC/jT c = 1.43, 
we obtain 7 = 9.6 mJ/mol- K 2 , small compared to 
that of non-magnetic counterpart, i.e. 18 mJ/mol-K 2 
for Lu(Y)Ni2B2C, possibly indicating a lower density of 
states, N(Ep), in H0NI2B2C. 19 A spectroscopic study, 
however, found that the density of states hardly changes 
within the borocarbide series (RNi2B2C). 20 Recently, El- 
Hagary et al. found a common correlation between the 
specific heat jump AC and the transition temperature 
T c among magnetic borocarbide superconductors, i.e., 
AC oc T 2 ^ which indicates the importance of the mag- 
netic pair breaking. According to Abrikosov-Gor'kov the- 
ory (AG), the exchange interaction between electrons and 
magnetic impurity atoms leads to nonconservation of the 
electron spin, affecting the formation of Cooper pairs. 
Assuming 7 ~ 18 mJ/mol-K 2 (same as that of Y or Lu 
based borocarbide) 2 ^, the ratio AC/jT c becomes 0.76. 
If we assume a 50 % suppression of T c from the nonmag- 
netic value of 16 K, the corresponding AG prediction is 
AC/ 7 T C ~ 1. 

It can be speculated that the origin of the small specific 
heat discontinuity results from anisotropic superconduc- 
tivity or a multi-band superconductivity as is the case 
for the two-gap superconductor MgB 2 < 18 i 23 ' 24 i 25 ' 26 In- 
deed, there is compelling evidence that the non-magnetic 
members of the borocarbides, Lu(Y)Ni2B2C, are highly 
anisotropic or probably nodal superconductors where 
there exist gap zeros on the Fermi surfaceiSL 2 ^^ The 
positive curvature in the upper critical field H C 2 of 
Lu(Y)Ni2B2C close to T c and the temperature depen- 
dence of H C 2 were successfully explained by using an ef- 
fective two-band models However, the gap anisotropy or 
multi-band feature is probably irrelevant to the anoma- 
lous value of AC/jT c in H0M2B2C because the reported 
thermodynamic ratio (= 2.3) of LuNi 2 B 2 G2I is much 
larger than the weak-coupling BCS value (=1.43) as well 
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FIG. 2: The top panel describes the specific heat of sample A 
as a function of temperature in several constant fields along 
[100] and the bottom panel represents the specific heat of 
sample B in fields along [110]. Both plots are on a semi- log 
scale 



as that of HoNi 2 B 2 C (= 0.76). 



III. MAGNETIC PHASE TRANSITIONS 
A. Magnetic transition at 5.2 K (Tjv) 

The top panel of Fig. 2 shows a semi-log plot of the spe- 
cific heat of H0NI2B2C (sample A) as a function of tem- 
perature at several magnetic fields B | [100]. The bottom 
panel of Fig. 2 describes the specific heat of sample B for 
B || [110]. The c-axis commensurate antiferromagnetic 
(AF) transition, labeled TV, is lowered with increasing 
field for both field directions. When field is higher than 
0.4 T, the AF peak becomes severely broadened, making 
the data difficult to interpret. The top panel of Fig. 3 
shows the B — T phase diagram for this transition to a 
Neel state. The critical temperatures T/v from C vs. T 
in constant fields (squares) and the critical fields Bm 
from the isothermal C vs. B at constant temperatures 
(circles; see Figs. 4 and 5) were plotted for both field di- 
rections. There is negligible in-plane Bm anisotropy be- 
tween [100] and [110] directions. The dashed line is the 
least-squares fit of B^ — A(l — T/Tjy) - 5 in the vicin- 
ity of Tjv, showing Ising-like behavior Experimental 
data were explained well with A = 1.41 ± 0.01 T. 

The bottom panel of Fig. 3 shows the hysteresis of 
the specific heat of sample B as a function of tempera- 



FIG. 3: The top panel describes Bm — T phase diagram for 
B || [100] (open symbols) and B || [110] (solid symbols). 
The circles represent the data from isothermal C vs. B and 
the squares, from C vs T in several fields. The dashed line 
describes the least square fit of A(l - T/T*) ' 5 with A = 
1.41 ± 0.01 T. The bottom panel shows the specific heat of 
sample B as a function of temperature in zero field. The 
circles represent the data with increasing temperature and the 
crosses, with decreasing temperature. Inset: the temperature 
offset of the sample as a function of temperature. 



ture. The circles describe the data with increasing tem- 
perature while the crosses, with decreasing temperature. 
There is no hysteresis in either the 6 K transition Tm nor 
the 5.5 K transition T*. In the AF transition, however, 
there is a clear hysteresis where the AF peak was moved 
from 5.21 K with increasing temperature to 5.16 K with 
decreasing temperature. The inset of the bottom panel 
shows the dc temperature offset of the sample as a func- 
tion of temperature. Depending on the direction of the 
temperature sweep, there occurs a sudden jump (or drop) 
in the dc offset at T)y, a consequence of the latent heat as- 
sociated with the transition. The rectangle-shaped hys- 
teresis adds additional strong evidence that the Tjv tran- 
sition is first order. This finding is consistent with the 
observation of magnetoelastic effects in H0M2B2C, where 
the length of the unit cell in [110] direction is shortened 
by about 0.19 % compared to its length in [110] direction 
below T/vf2 3 Despite the square-root behavior evident in 
Fig. 3, the first-order nature of the magnetic AF transi- 
tion disagrees with the Ising models We note that the 
peak intensity is smaller for the specific heat with de- 
creasing temperature. That could be an artifact from 
ac calorimetry because the temperature response of the 
sample is not an ideal triangular shape, rather a distorted 



4 



B || [100] 




B || [110] 



FIG. 4: Isothermal specific heat of sample B as a func- 
tion of magnetic field at several temperatures for B || [100]. 
Prominent peaks and shoulders were marked by arrows. The 
peaks and shoulders are more evident on an expanded scale. 
The transition from superconducting state to normal state is 
marked as T c . The dotted lines are guides to the eye. 
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FIG. 5: Isothermal specific heat of sample B as a function of 
magnetic field at several temperatures for B || [110]. Promi- 
nent peaks and shoulders were marked by arrows. The peaks 
(or shoulders)at each temperature were determined on an ex- 
panded scale. The transition from superconducting state to 
normal state is marked as T c . The dotted lines are guides to 
the eye. 



one when it goes through a phase transition with latent 
heat M. 



B. The 5.5 K magnetic transition (T*) 



The 5.5 K magnetic transition (T*) is as sharp as the 
AF transition, for sample A at least, indicating it could 
also be a first order phase transition (see Fig. 2). The 
magnetic field dependence of the specific heat is also sim- 
ilar to the TV counterpart, i.e. the critical temperature 
decreases with increasing magnetic field. Initially, the 
peak intensity at T* becomes stronger with increasing 
field while that of the AF transition decreases monoton- 
ically, transferring some of its entropy to the T* transi- 
tion. Fig. 6 shows the B* ~ T phase diagram both for 
B || [100] (open) and B || [110] (solid), where B* is the 
critical field corresponding to the T* transition. The in- 
plane B* anisotropy between [100] and [110] directions is 
small, similar to the Sjv transition. The B* temperature 
dependence near T* was explored in terms of the anti- 
ferromagnetic theory by Fisher ^ as was done in the TV 
transition. The suppression of the critical temperature 
was explained reasonably well by B* = A(l - r/T*) - 5 
with A = 1.39 ± 0.01 T (dashed line), a value similar 
to that of the Tm analysis. 



What could be responsible for the 5.5 K transition? 
Neutron scattering showed that all three magnetic struc- 
tures, i.e. c-axis commensurate AF magnetic structure, 
c-axis AF helical structure, and a-axis incommensurate 
structure, coexist between 5 K and 6 In the pre- 

vious section, we argued that the 5.2 K transition is due 
to the magnetic transition to the c-axis commensurate 
AF structure (TV) and is a first order phase transition. 
All the similarities between the 5.5 K transition and the 
5.2 K transition point toward associating T* with the c- 
axis incommensurate AF structure. Then, the next ques- 
tion is "Is it a first order phase transition?" The specific 
heat at zero field did not show any noticeable hysteresis 
with temperature at T* (see Fig. 3). In the dc tempera- 
ture offset (see the inset of Fig. 3), however, there may be 
a small hysteresis at T* , but the feature is well within the 
scattered data. Coexistence of the phases also suggests 
that it is 1st order. At this point, it is not clear if the 
5.5 K magnetic transition is a first order phase transition 
or not. 
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FIG. 6: B* — T phase diagram for both field directions: open 
symbols describe the data for B || [100] and solid symbols, for 
B || [100]. The squares represent the data obtained from C 
vs. T while the circles, from isothermal C vs. B data. The 
dashed line describes the least square fit of A(l — jyy*) - 5 
with A — 1.39 ± 0.01 T. 
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C. The 6 K magnetic transition (Tm) 

The 6 K magnetic transition (Tm ) is known to be the 
onset of a long-range magnetic order. The precise nature 
of the Tm transition is not well characterized, although 
recent work2& argues that it represents the onset of the a* 
modulated phase. In this section, we present the specific 
heat as a function of temperature both for B || [100] and 
B || [110] directions and confirm that the 6 K transition is 
due to the onset of the a-axis incommensurate magnetic 
structure. 

Fig. 2 shows a dramatic difference in the critical tem- 
perature Tm with magnetic field directions. For B | 
[110], the critical temperature scarcely changes with in- 
creasing magnetic field below 0.5 T in agreement with 
Detlefs et ali^ii and Du Mar, et aliS. In higher fields, it 
decreases rapidly as is expected for an antiferromagnetic 
transition and the transition shape becomes broadened, 
making it hard to interpret. For B || [100], the criti- 
cal temperature increases with increasing magnetic field- 
while keeping its transition shape. The anomalous in- 
plane anisotropy of Tm strongly suggests that the ori- 
gin of this transition is very different from the other two 
magnetic transitions (TV, T*) and is, therefore, due to 
the onset of the a— axis incommensurate magnetic struc- 
ture. Detlefs et ali2& recently hypothesized that the a* 
magnetic structure is related to the Fermi surface (FS) 
nesting along [100] direction. Since a distortion along 
[110] direction is more disruptive to the nesting feature 
than a distortion along [100] direction, it results in the 
absence of the a* phase for H0M2B2C and DyN^E^C 
at low temperature and low magnetic field where local 
magnetic moments are aligned along [110] directions. If 
the a* phase is the ground state for 5.5 K < T <6.0 K in 
H0M2B2C, a magnetic field applied along [110] will dis- 



FIG. 7: The top panel blows up the specific heat of sample A 
vs. temperature near Tm in several magnetic fields when B | 
[100]: from left to right, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 
1.0, 1.5 T. The bottom panel shows a scaling behavior of T* 
transition where x-axis is the reduced temperature 1 — T /Tm 
and y-axis is the normalized specific heat C /Ct m where Ct m 
is the specific heat at Tm- For clarity, selective data were 
shown. 



rupt the magnetic phase while the field along [100] will 
not. The anomalous magnetic field dependence of specific 
heat is consistent with this scenario. The enhancement 
of the a* phase for B || [100], however, is beyond this 
explanation. 

The top panel of Fig. 7 shows the specific heat data 
around Tm in several constant magnetic fields from T 
to 1.5 T for B || [100]. For clarity, selective data are 
shown in the bottom panel, where x-axis is the reduced 
temperature (1—T/Tm) and y-axis is the specific heat di- 
vided by the specific heat at Tm, i.e. C/Ct m - All of them 
collapse onto each other, showing a scaling behavior. The 
specific heat near a critical point diverges logarithmi- 
cally in the two-dimensional Ising model while it diverges 
more strongly, as a power law in the three-dimensional 
Ising models Fig. 8 describes the zero-field magnetic 
specific heat as a function of the reduced temperature 
t(= 1 — T/Tm) on a semi-log scale. The magnetic specific 
heat was obtained by subtracting the lattice and elec- 
tronic contributions: Cm = C — jT — (3T 3 where we used 
7 and values of TmN^E^C,? 8 The higher temperature 
side of Tm was analyzed in terms of 2D and 3D Ising 
model. The dashed line is the best fit of 3D Ising model 
with a functional form of A | -r | 0-1 + B. The solid line 
is from 2D Ising model of — C log\r\ + D. Both 2D and 
3D Ising models can explain the data over two decades 
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FIG. 8: Magnetic specific heat of sample A as a function of 
the reduced temperature r(= 1 — T/Tm) in zero field on a 
semi-log scale. The dashed line is the least square fit of 3D 
Ising model, i.e. A |r|~ 01 + B with A = 9.92 J/mol K and 
B — —6.68 J/mol K. The solid line represents the least square 
fit of 2D Ising model, i.e. -C log\r\+D with C = 3.0 J/mol K 
and D = 2.9 J/mol K. 
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of temperature range, i.e. 3 x 1CP 3 < r < 3 x 1CP 1 . 
When temperature is close enough to Tm(— 5.98) K, 
however, the specific heat data can be explained better 
with the logarithmic function than the power-law depen- 
dence. The logarithmic singularity may be interpreted 
as a manifestation of two-dimensional spin structure in 
H0M2B2C where the Ho 3+ local moment is confined to 
the Ho-C basal plane for temperatures below 100 KjS We 
note that the value of the coefficient of the logarithmic 
term, C = 3 J/mol- K, is of the order of magnitude found 
in the exact theory of 2D Ising antiferromagnets^ The 
lower temperature side of Tm is severely contaminated by 
the adjacent magnetic transition T* , rendering analysis 
difficult. 



D. Magnetic phase diagram 

It has long been known that there is an extreme mag- 
netic anisotropy associated with the crystalline electric 
field (CEF) splitting of the Hund's rule ground state 
for the magnetic members of the borocarbide family 
RNi 2 B 2 C (R=Er, Tb, Ho, Dy)iLia The R 3 + local mo- 
ment is confined to the R-C basal plane for temperatures 
below roughly 100 K, i.e. temperatures well above the 
magnetic ordering temperatures. In addition, a strong 
in-plane anisotropy has been observed, leading to the lo- 
cal moments essentially being confined to either the [100] 
(R = Er, Tb) or [110] (R = Ho, Dy).4i Therefore, it be- 
came necessary to study physical properties as a function 
of magnetic field direction. 

Fig. 9 summarizes the magnetic phase diagram of 
H0M2B2C both for B || [110] (upper panel) and for 
B || [100] (lower panel) directions based on specific 



FIG. 9: B — T magnetic phase diagram both for B || [110] 
(upper panel) and for B | [100] (lower panel ) directions. 
The triangles describe the c— axis commensurate AF transi- 
tion, the circles the c— axis incommensurate AF transition, 
and the squares the a— axis incommensurate transition. The 
crosses in the upper panel represents an additional phase line 
for B || [110], most readily seen as a shoulder in the 4.8 K 
and 5.0 K in Fig. 5. Dashed lines are guides to the eye. The 
insert in the bottom panel is the full Bm —T phase diagram 
up to 6 T for B || [100]. 



heat as a function of temperature and magnetic field. 
The 5.2 K transition, TV, is due to the commensurate 
c— axis AF structure and found to be a first order phase 
transition. The critical field Bm (triangles) shows lit- 
tle in-plane anisotropy. The 5.5 K transition, T*, is 
due to the c-axis oscillatory magnetic structure. The in- 
plane anisotropy of B* (circles) is small and the nature 
of the transition is still controversial. The 6 K transi- 
tion, Tm, is ascribed to the onset of the a-axis oscilla- 
tory magnetic structure. The in-plane anisotropy of Bm 
(squares) is very surprising: the critical temperature 
decreases with increasing magnetic field for B || [110] 
while it increases, for B \\ [100]. The critical tempera- 
ture Tm initially increases with magnetic field at a rate 
of 0.48 K/T, then it starts to bend around at about 1.5 T 
and decreases (see the inset of the bottom panel of Fig. 9). 
The decrease of the critical temperature is expected for 
antiferromagnet. 31 The increase, however, is difficult to 
understand. 

The magnetic phase diagram built from specific heat 
data is consistent with the viewpoint of three distinct 
magnetic transitions between 5K < T < 6K at zero 
fieldAii At low fields, the critical temperature Tm 
scarcely changes with magnetic field for B \\ [110] in 



agreement with Du Mar et al. 9 while the increase of Tm 
for B || [100] is more evident, probably due to high- 
quality sigle crystal used in this study. In the high- 
field regime, our data show a dramatic difference in Tm 
between the two field directions, i.e., B || [100] and 
B || [110] as is evident in Fig. 9 while previous studies 
drew a conclusion of similar behavior between the two 
directions. For B \\ [110], Ref 9 shows only two transi- 
tions below 5 K, i.e., T/v and T*. In contrast, our data 
clearly show three transitions at the same temperature 
range. A close examination reveals that the T* marked 
transition below 5 K in Ref 9 corresponds to Tm in our 
phase diagram, indicating that the real T* transition is 
missing in that phase diagram. 

There has been a speculation of a fourth magnetic 
phase based on resistance measurements where there oc- 
curs a slope change at 3.8 K for B || [100] In our bulk 
measurement, the additional feature was not observed, 
suggesting that the 3.8 K feature is an extrinsic property 
sensitive to the surface state. Instead, we found an ad- 
ditional phase line for B || [110], most readily seen as a 
shoulder in the 4.8 K and 5.0 K data in Fig. 5. It sepa- 
rates the Bm and the B* lines and is merged to Bm at 
4 K (crosses in the upper panel of Fig. 9). The origin of 
this feature has yet to be elucidated. 



IV. MAGNETIC FIELD-ANGLE HEAT 
CAPACITY 
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FIG. 10: Heat capacity of sample B as a function of field angle 
at several temperatures. The magnetic field is rotated within 
the basal plane and is measured with respect to a-axis. The 
detailed condition is marked in each responsible panel. 



In non-magnetic borocarbides, the low temperature 
heat capacity directly measures the electronic density 
of states^ In magnetic systems, however, the magnetic 
specific heat dominates. For example, the electronic and 
lattice components of the specific heat account for less 
than 5% of the total below 8 K in HoNi 2 B 2 C. Conse- 
quently, in HoNi2B 2 C, the field-angle dependent heat 
capacity used to investigate the superconducting gap na- 
ture in the non-magnetic systems, here mainly explores 
the magnetic structure. In this section, we present the 
in-plane field-angle heat capacity of the Ho-based boro- 
carbide and discuss the data based on the phase diagram 
built in the previous sections. 

Fig. 10 shows the low-field heat capacity vs field angle 
measured with respect to the a— axis of H0NI2B2C at 
several different temperatures. The periodicity of the 
peaks is 90 degrees at all temperatures, indicating that 
the peak in the specific heat is simply due to the in-plane 
anisotropy between [100] and [110] directions. Canfield 
et al. observed magnetization modulation as a function 
of magnetic field angle at 2 K4i The oscillation feature 
was then interpreted in terms of metamagnetic states, 
i.e., (U) for B < B c i, (TTI) for B cl < B < B c2 , (tH) 
for B c2 < B < B c3 and (ftt) for B > B c3 . Here the 
arrow f is a moment along the [110] axis, J, is a moment 
along the [110], and — > is a moment along the [110] axis. 
In our magnetic phase diagram (Fig. 9), B c \ corresponds 
to Bn, B C 2 to B* , and B c3 to Bm respectively. 



At 2.5 and 6.1 K, the low-field heat capacity has max- 
ima for the field along < 100 > directions. The area un- 
der the field-angle heat capacity is proportional to mag- 
netic entropy change, indicating that there is more mag- 
netic disorder for the field along < 100 > than along 
< 110 >. The field-angle entropy modulation can be ex- 
plained by the fact that the net moment of Ho 3+ ions 
is directed along < 110 > directions in the commensu- 
rate antiferromagnetic phase. 35 At 4.8, 5.2, and 5.6 K 
where helical magnetic phases appear, the peak positions 
shifted by 45° to < 110 > directions, indicating that the 
magnetic moments are preferably aligned along < 100 > 
directions. The peak intensities are also asymmetric with 
field angle, which may be related to the oscillatory mag- 
netic structures observed 3 *^ in this temperature range. 

Fig. 11 shows the field-angle heat capacity at 6.1 K in 
0.5, 1, 2, and 4 T. According to the magnetic phase dia- 
gram of H0M2B2C (see Fig. 9), the heat capacity should 
exhibit monotonic 4-fold angular oscillation at high fields 
because there is no transition other than the Tm tran- 
sition nearby this temperature. The peaks along [100] 
at 0.5 T, however, were split into two peaks at 1 T and 
the minima along [110] were totally flattened out. At 
2 T, the distance between the two split peaks becomes 
narrower and the flat minima become broad maxima. In 
4 T, the split peaks merge and show a delta function like 
peak and the broad maxima at [110] in 2 T returns to a 
broad minima as was in 0.5 T. Two satellite peaks appear 
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FIG. 11: Heat capacity of sample B as a function of field 
angle at 6.1 K in 0.5, 1, 2, 4 T. The circles describe the data 
with rotating magnetic field clockwise and the crosses, with 
rotating field counter clockwise. 

at ±12° of [100] peaks. 

The seemingly anomalous peak splittings along [100] 
may be explained as following. In 1 T, Tm moves from 
3.7 K for B || [110] to 6.5 K for B || [100] or T M = 
3.7 + 0.031a where a is the in-plane field angle measured 
against the a-axis. To move the Tm through 6.1 K means 
a field angle of 0.4 K/0.031 K-deg" 1 w 12.5°, which is 
consistent with the peak positions at 10± 1.5. At 2 T, the 
T M moves from K for B || [110] to 6.57 K for B || [100], 
predicting 6.3° of peak splitting. Experimentally, the 
observed splitting was 6 ± 1.5°. The larger the field, the 
faster the peak has to move in angle. Even though the 
anomalous peak splitting could be accounted for in terms 
of a simple linear relation between Tm and field angle a, 
the interchange between minima and maxima at [110] 
and the satellite peaks in 4 T have yet to be understood, 
suggesting that the magnetic phase in high fields is not 



as simple as was originally envisaged. 

In Fig. 11, the circles describe the data when the 
in-plane magnetic field was rotated clockwise and the 
crosses, when the field was rotated counterclockwise. In 
this measurement, the angle spacing between two data 
points is 3°. At 0.5 T, the two data sets do not show hys- 
teresis in phase, but show hysteresis in amplitude with 
field direction. Above 1 T, there is a systematic shift in 
phase by ~ 3° with field direction, indicating that the 
related transition may be a first order magnetic transi- 
tion. The experimental error bar in determining the field 
angle is ±1°, which is less than that of the angle shift 
in peak positions. The hysteresis with field- angle direc- 
tion may be another manifestation of the first order na- 
ture of the metamagnetic transition from (TT~~ to (TTt) 
phases, which were observed in neutron-diffraction study 
by Campbell et al42 



V. SUMMARY 

We have studied the magnetic superconductor 
H0N12B2C via specific heat as a function of tempera- 
ture, magnetic field, and magnetic field angle. The small 
value of the specific heat discontinuity at T c indicates 
that Ho 3+ ions act as pair breakers as was predicted in 
Abrikosov-Gor'kov theory. The temperature hysteresis 
of the zero-field specific heat at Tjy provides direct evi- 
dence that the magnetic AF transition is of first order, 
accompanying a change from tetragonal to orthorhombic 
structure. The similarity of T* transition to the Tjv tran- 
sition indicates that it is due to the c* spiral structure. 
The anomalous in-plane anisotropy of Tm transition was 
explained in terms of the 0.6 a* Fermi surface nesting 
feature. 
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